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Abstract 

D , 

. An effective theory to treat the dense nuclear medium by the perturbative 

expansion method is proposed as a natural extension of the Heavy Baryon 
^ ■ Chiral Perturbation Theory (HBChPT). Treating the Fermi momentum scale 

\ as a separate scale of the system, we get an improved convergence and the 

conceptually clear interpretation. We compute the pion decay constant and 
the pion velocity in the nuclear medium, and find their characters different 
from what the usual HBChPT predicts. We also obtain the Debye screening 
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scale at the normal nuclear matter density, and the damping scale of the pion 
wave. Those results indicate that the present theory, albeit its improvement 
over the HBChPT, has the limitation yet to go over to the medium of about 
1.3 times of normal matter density due to the absence of the intrinsic density 
dependence of the coupling constants. We discuss how we overcome this 
limitation in terms of the renormalization method. 

1 Introduction 

Low-energy nuclear processes that involve a few nucleons can be described in a sys- 
tematic and consistent way, thanks to the heavy-baryon chiral perturbation theory 
(HBChPT) which has been initiated by the pioneering work of Weinberg . This 
noble approach has been apphed to many low-energy processes with great successes; 
see, for example, Refs. |2l El IH E] ■ 

Motivated with these successes, many already applied the chiral perturbation 
theory to dilute systems [3 13 IHl El UH] and also to systems around the saturation 
density When the density and the corresponding Fermi momentum pp 

increases, however, the convergence of free-space HBChPT becomes questionable. 
This is because HBChPT is a derivative expansion scheme with the expansion pa- 
rameter Q/A^; Q stands for the typical momentum scale (and/or the pion mass) of 
the process, and ~ 4-7? /,r — 1 GeV is the chiral scale. In nuclear matter, Q ^ pp 
even for small fluctuations whose typical energy-momentum scale is much smaller 
than Pp. For such fluctuations, the contribution of the nucleons whose momenta 
far from the Fermi surface is however irrelevant. It is thus desirable to build an 
effective field theory (EFT) by integrating out such ^^massive degrees of freedom", 
\p\—pp > A, where A is the cutoff of the theory. The resulting EFT has two expan- 
sion parameters, Qext/-^x Pf/-^x^ where Qext represents the energy-momentum 
scale of the external probes and/or the pion mass, Qext < A. The resulting EFT is 
thus meaningful only when Qext < A < pp and pp < A^. These expansion parame- 
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ters will emerge naturally in the power counting rule of the matrix elements of the 
physical processes, and the perturbative expansion will follow this counting rule as 
is the case in the HBChPT where only the expansion scale is Qext < A^. As we shall 
see later, the new expansion scale pp < could be threatened to approach to one 
when the density of the nuclear medium becomes large. But as we shall show later, 
the power counting rule shows that the matrix elements depends on the pp/A^ to 
twice of the number of loops that consist only of nucleon lines. This guarantees 
the rapid convergence in the dilute situation. In this work, however, we limit our- 
selves to a modest case where Qext ^ ^ A^, which allows the expansion both 
in Qext/ A^ and in ppj A^. This is the region that can be treated by the HBChPT. 
However, the EFT that we are developing here has two important advantages over 
the HBChPT. Firstly, the EFT has only relevant degrees of freedom (integrating out 
the irrelevant ones), which is the key aspect of EFTs that enables us to describe the 
nature systematically and consistently in a transparent manner. Secondly, since the 
convergence with respect to powers of Qext/ A^ is much faster than that of pp/A^, 
substantial reduction of the relevant diagrams can be achieved. 

The coefficients of the resulting EFT Lagrangian, i.e., low-energy constants 
(LECs), will be in general different from those in the vacuum and must be de- 
termined from QCD. At this moment, we do not have any suitable ways to match 
our EFT directly with QCD. Instead we fix them by matching the EFT with the 
HBChPT; the Fermi momentum of the matching point should be upper-bounded 
by the validity region of HBChPT and lower-bounded by that of our EFT. Once 
the LECs are fixed, we determine the pF-dependences of the physical quantities by 
including the (nucleon-)loop corrections. 

Before going further, we acknowledge that the idea of building EFT near the 
Fermi surface is not new at all. For example, the familiar Landau-Migdal Fermi 
theory is one of such. The work of Schwenk et al. also is worthy to be noted, 
where the notion of the Viow-k [El has been applied to derive the parameters of 
the Landau-Migdal theory. Especially our work shares much in common with the 
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high density effective theory (HDET) developed by Hong J3] and by Schafer [TE] . 
which is an EFT for quark matter at high-dense region where the matching with the 
(perturbative) full QCD is possible. Compared to the HDET, our prime interests lie 
in the intermediate density (rather than the asymptotic) region where the (almost 
massless) Goldstone bosons exist and the relevant Fermion degrees of freedom are 
still nucleons. However, our approach may extend its applicability to the region 
where the chiral quarks are the relevant degrees of freedom. Our present study 
might be relevant for the recent pionic atom experiment at GSI with heavy nuclei, 
where the pion decay constant is interpreted to decrease that may signal the partial 
chiral symmetry restoration in medium |17j . 

We present in this paper our first step to make the in- medium EFT. Our explicit 
degrees of freedom are nucleons near the Fermi surface, |p| — < A and the 
Goldstone bosons. As is claimed recently by one of us (MH) ^TF , the role played by 
vector mesons at high density should not be underestimated to realize the important 
phase of matter at the transition, the Vector Manifestation (VM) JHH^Ol- Especially, 
we believe it is very important to understand how this phase is connected from the 
normal nuclear matter. Our model indeed stems from the purpose of its realization. 
In the present report, we are however much modest to consider only the pion degrees 
of freedom. The underlying rationale lies in the decimation process at hand. Our 
theory will allow the shallow excitation nearby the Fermi surface up to, say, the 
pionic mass scale. Those modes beyond it will be integrated out from the beginning. 
In a sense, our model will meet the limit of the applicability if the actual mass of 
certain integrated-out degrees of freedom becomes lighter as the density increases. 
This phenomena will emerge in our EFT when the parameters do not satisfy the 
naturalness. In this respect, it will be sought in the forthcoming work to include the 
chiral partner of pion in the dense limit. 

It is worth mentioning some of our ansatze. We are dealing with a dilute system 
in which loop corrections to the 4-Fermi type interaction are suppressed by some 
powers of pp/A^. Furthermore, the most important channel is assumed to be the 
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forward scattering one as in the Fermi liquid pi], which builds Landau's Fermi liquid 
theory [22] • In the present analysis, we neglect BCS phenomena. The incorporation 
of the BCS channel is postponed to the next step. 

This paper is organized as follows: We give a brief illustration to define the field 
expressing the fluctuation mode of the nucleon in section |21 and list first few terms in 
section 121 In section |31 we build the power counting rule for our effective Lagrangian 
and list more higher order terms. The in-medium modification of the vector corre- 
lator is studied shortly in section El In section (HI we study the corrections to the 
pion decay constants and the pion velocity using our EFT. The bare parameters for 
the decay constants will be determined through the matching between our EFT and 
the HBChPT. In section [71 we consider the mass(-like) terms and check the change 
in quark condensate due to mass term in our theory. We summarize our works in 
section jSl 

2 Nucleon and its fluctuation flelds 

In this section, we will define the field corresponding to soft modes of the fluctuation 
around the Fermi surface and obtain the higher order correction terms in the present 
model scheme. 

Let us begin with the free Lagrangian for the nucleon without interaction at 
non-zero density: 



where /z is the baryon chemical potential and rriN the nucleon mass. Fermi momen- 
tum is related to the mass mi^ and the chemical potential /i as in 



In the dense medium of interacting nucleons, we should use the in-medium nucleon 



C = i) {id^-i^ -mN + /x7o) , 



(2.1) 




(2.2) 
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mass, say m*^, instead of the one in vacuum But for the notational simphcity, 
we refer to as the in- medium mass m*^ throughout this paper. 

In order to consider the fluctuation mode around the Fermi surface, we introduce 
the following field ipi. 



(2.3) 



Vp 



where the Fermi velocity vp is related to the Fermi momentum pp as 

1 ^ 



vf = -Pf 
/i 



(2.4) 



Simply the sum of vj means the angle sum on the Fermi surface. By using the field 
iPi{x;vf) in Eq. ()2.3j) . the Lagrangian in Eq. 1)2.11) is rewritten as 



= XI ^1 (^'9/,7^ - Pf ■ 1 - rriN + yU7o) 'ipi 



(2.5) 



Vp 



where 7 = (7^, 7^, 7^). To delete off-diagonal components of the 7 matrices between 
'01 and ipi we introduce the following unitary transformation: 



ipi{x; vf) = Uf{vf) ip2{x] vf) 



(2.6) 



where 



Uf[Vf) = cosfc* z — smt^ = cos 6^ smb* 



Vf 



Pf 



with Vf = \vf\ and 9 satisfying 



cos = \ — , sm (J 



I n — mjv 



(2.7) 



f2.8) 



' 2/i V 2/i 

By using the field 4'2{x; Vf) in Eq. ()2.(ij) . the Lagrangian in Eq. ()2.5|1 is rewritten as 



Vp 



i9o7° + idjV^p + idj ( 7^ + \ + /"To - 



'2 , 



(2.9) 
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where summations over j are implicitly taken and Vp denote components of vp as 
Vp = {vF,Vp,Vp). From now on, the indices i and j run over 1, 2 and 3. Two 
7 matrices representing the parallel and the orthogonal components to the Fermi 
velocity, 7| and 7j_, are defined by 

7fi = i<(^F-7) , (2.10) 

Vp 

li = l' -^vUvF-l) . (2.11) 

Vp 

Now, we decompose the '?/'2 field into two parts as 

tP±{x;vp)=P±Mx-,VF) , (2.12) 

where the projection operators P± are defined as 

P, . . (2.13) 

and ip- represent particle and antiparticle fields in low energy, respectively. By 
using the fields 4'±{x; vp), the Lagrangian in Eq. ()2.9p is rewritten as 



C = J2U+ vp id^ ^+ + ^. {-vp id^ - 2/i) v^. 



Vp 



(2.14) 



where 



VP = {l,vp), vp = {l,-vp). (2.15) 



The Lagrangian ()2.14|) shows that the field expresses the fiuctuation mode of the 
nucleon around the Fermi surface, and the field ip- corresponds to the anti-nucleon 
field carrying the effective mass of 2/i which is large compared with the momentum 
scale we are considering. Then, we can solve out the ■?/'_ field in studying the low- 
momentum region to obtain the Lagrangian including the field. It is to be noted 
that doipj^ is of same scale as diip+, thanks to the projection. 
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In the next section, we will construct the general Lagrangian in terms of the V'+ 
field. Here to have an idea on the higher order terms, we solve out the ip- field in 
the Lagrangian ()2.14j) at tree level. Using equation of motion for the ip- field 



7l + ^7f 



(2.16) 



we obtain 



Vf 



^/^ n=0 



2/i ^ 



(2.17) 



with 7 = 7± + ^^7||- This shows that the terms in the right-hand-side (RHS) other 
than the first term are suppressed by Q/ii, where Q denotes the typical momentum 
scale we are considering. Therefore, the first term is the leading order term, which 
leads to the following propagator for the fiuctuation field '0+ 

1 



iSpip) 



-Vf ■ p - iepo 



(2.18) 



where e — > +0 and P+ is the projection operator defined in Eq. ()2.13p . Other terms 
are regarded as the higher order terms. Thus, the leading correction to the first 
term of the RHS is already suppressed considerably at any density. 



3 Effective Lagrangian 

In this section, we will present our effective Lagrangian. As stressed in Introduction, 
we are considering the moderate density region. So the mesons other than the pion 
are expected to be still heavier than the pion, i.e., the pion is the only light mesonic 
degree of freedom realized as the pseudo Nambu-Goldstone (NG) boson associated 
with the spontaneous chiral symmetry breaking. In the present analysis, therefore, 
we build an effective Lagrangian including the field expressing the fiuctuation of the 
nucleon around the Fermi surface, together with the NG boson field, based on 
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the non-linear realization of chiral symmetry by assuming that all the other mesons 
are integrated out. 

As we illustrated in Introduction, we are considering the nucleon carrying the 
momentum inside the thin shell of size A|| around the Fermi surface. Then, the 
momentum of the fluctuation parallel to the Fermi velocity at the reference point 
is restricted to be smaller than A||. Furthermore, we only consider the situation 
where the momentum of the fluctuation mode perpendicular to the Fermi velocity 
is also smaller than A_l ~ A|| by assuming that effects from mode carrying large 
perpendicular momenta are already integrated out and included in the parameters 
of the effective Lagrangian. This idea is based on Landau's Fermi liquid theory 
where the forward scattering channel is most important PT] and which works well in 
nuclear physics j23- And such a scheme is being used in HDET for quark matter [T2j. 
Thus, in the present analysis, we divide the thin shell around the Fermi surface into 
small boxes, each of which is characterized by the point on the Fermi surface with 
the characteristic vector, vp. Since the size of each box is A||A^, the momentum of 
pion is restricted to be small and the nucleon interacting with pion cannot escape 
from the box. 

Terms without fermion 

The pion is introduced as the NG boson associated with the chiral symmetry 
breaking of SU(2)l x SU(2)r SU(2)v through the matrix valued variable ^ 



where vr = I]a=i ^aT'a denotes the pion fields, and the parameter corresponding to 
the temporal component of the pion decay constant. Here the wave function renor- 
malization of the pion field is expressed by ^ |21] . Under the chiral symmetry, 
this ^ transforms nonlinearly as 



(3.1) 



(3.2) 
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where h(7i , gn, gi^) is an element of SU(2)v to be uniquely determined. There are 
two 1-forms constructed from ^: 

1 



a!^^-(l)^^^t_2)M^t.^j ^ (3.3) 
^ 1 {V^^ . e + V^^^ ■ , (3.4) 

where V^^ and V^j,^"^ are defined as 

V,^^ = d,e + z^^C^ , (3.5) 

with £ and TZ being the external gauge fields for the chiral SU(2)l x SU(2)r sym- 
metry. These 1-forms transform as 

"a 9r, 9l) a a h\T^, 9r, gh) , 

h{Tr, gK, gh) h\T^, Qr, 9l) - -h{TC, g^, g^) d^" h\-n, g^, gh) ■ (3.6) 

By using these 1-forms the general form of the chiral Lagrangian including pion 
fields can be constructed. The Lagrangian at the leading order of the derivative 
expansion is given by 

, 2 



AO 



(f*) u,u. + (f^F:) (g, 



where = (1,0) indicates the rest frame of medium, and F^ and F^ denote the 
parameters corresponding to the temporal and spatial pion decay constants. In our 
notation {0\dnA^\7r) oc F^Pq — F^p^ = for on-shell pions by axial vector current 
conservation [23] so the pion velocity is expressed as = F^/F^. It should be 
noted that since F^ is the wave function renormalization constant of the pion field, 
the chiral symmetry breaking scale is characterized by 

Ax ~ 47rF*. (3.8) 
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Fermion kinetic term 

Now, let us consider the terms of the effective Lagrangian including the field for 
expressing the fluctuation mode of the nucleon around the Fermi surface. In the 
following, for notational simplicity, we use \1/ for expressing the fluctuation mode 
which was expressed by ip^ field in the previous section. Under the chiral symmetry 
the fluctuation mode \l/ transforms as 

^^/i(7r,(7R,(7L)^ , (3.9) 

where h{n,g^,gi^) G SU(2)v. Then the kinetic term of the \1/ field is expressed as 

C^in = T.^VPtD,^ , (3.10) 

Vp 

where the covariant derivative is defined by 

= - zav^) ^ , (3.11) 

with avfj, defined in Eq. ()3.4p . 



Interactions of fermions with mesons 

We consider the interactions of \1' and ^ to the vr field. The possible forms of 
the interaction among ^ and defined in Eq. ()3.3|) are given by 

^A = J2[^ f^AO^ivF ■ iHoiA^ + i KAW^li^f^aAi^ + i f^A±W±75aAi^] , (3.12) 

Vp 

where kao, i^a\\ and ka± are dimensionless real constants, and 7|j and Yj_ are defined 



\1/ are 



in Eqs. (ETT|l and TTT^ . *^ Note that -^a\^, -^v^aAi^ and ^) 

prohibited by the parity invariance. 

"^^It may be useful to compare the interaction terms in Eq. H3.12(l with those reduced from a 
simple interaction among pion and nucleons given by 
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Four Fermi interaction 

For constructing four- Fermi interactions, tlie following relations are convenient: 



^75^ = , ^7570^ = , 

'I' [7o , 7i] ^ = , ^ [7i , 7j] ^ = '^i^ijk ^157^"^ , 

[7. , 7 .] 75^ = , ^ [70 , 7i] 75^ = 2^'7i75^ • (3.13) 



Although the excitation momentum of the nucleon included in each pair of par- 
ticle and hole is restricted to be small in the present analysis, Fermi velocity of 
one pair can be different from that of another pair. Thus, using the relation in 
Eq. ()3.13|) . we can write the four-fermion interactions as 




^'l'7M75rvi>| (5]W^^7^^5fvi>| . (3.14) 



where r/2 is the generator of SU(2). It should be noticed that the above four- Fermi 
interactions allow two pairs of particle and hole, each of which resides in different 



By performing the same procedure as in section |21 the above term is reduced to 



{vF ■ 7)75aAo 



( i 






7l 




V M 







where 7^ and 7iJ are defined in Eqs. H2.1f |l and (|2.1U|) . Then, when we match the interaction terms 
in Eq. I|3.f 2|l with the above terms, we obtain the foUowing matching conditions: 



ruN 

KAO = HA\\ = 9A , KA± = gA 
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box. Here we introduce the dimensionless parameters by rescaling the dimensionful 
parameters by the temporal component of the bare pion decay constant. This is 
natural since we are constructing the EFT in the chiral broken phase |26j . 



4 Power counting 

In this section, we first present the power counting theorem in our EFT, which 
provides that the terms listed in the previous section are actually the ones of leading 
order. Then, we provide the next order terms according to the power counting. 

4.1 Power counting 

Let us consider the matrix element M with external nuclear lines, E^^ external 
TT lines and Ee lines corresponding to the external fields such as and IZ^ in 
Eq. (j3.5|) . Here, for simplicity we assume that all the mesonic external fields carry 
the dimension one and that all the nucleonic external fields the dimension 3/2. 

To derive a power counting rule for M, we first observe that the pion and nu- 
cleon propagators that carry the momenta of order of Q are of order of and 
respectively. Furthermore, we can count loops that consist at least one pion 
propagator as of order of / d!^q ~ vr^Q"^, where the factor vr^ is due to the angle 
integration. In fact, these observations can also be applied to free-space HBChPT. 
In our EFT, the integral is realized as 



since the summation over vp and the integral over the perpendicular space corre- 
sponds to the area of the Fermi surface. For the integrations over the momentum 
parallel to the Fermi momentum /|| and the energy Zq, we assign / dl^dlw ttQ"^, 
'^^We perform the Fourier transformation for each operator associated with the external field. 




(4.1) 
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where the overall factor vr can be understood either by the angle integration in 2-D 
space or by the residue calculation for the integration over Iq. 

Recalling the fact that each loop is accompanied by (27r)~^, we are thus led to 
the following equation for the counting of a Feynman diagram expressing M, 



M oc g-^^-^- ■ { S-^ i ?#1 Q' , (4.2) 



where {In) denotes the number of the pion (nucleon) propagators, L the number 
of loops, and d is the sum of the number of derivatives at z-th vertex di. Ln is the 
number of fermionic loops, each of which carries independent Fermi velocity vf, i.e., 
number of summations over vp. We note that we are working in the chiral limit, 
so that the pion mass does not appear in the above formula. Using the techniques 
developed in HBChPT, we can simplify the above equation as 

2L]v 



M^mV^^I . (4.3) 



with *^ 



E 



u = 2-(^ + Ee)+2{L-Ln) + Y1 



ft ' 

u, = di + ^ + ei-2 . (4.4) 

Here (cj) is the number of nucleons (external fields) at i-th vertex. For the 
detailed derivation of the above counting rule, we refer Ref. jIll2ZlEl- Note that we 
have replaced 4:71 by A^. 

It is worth while mentioning on the comparison of our counting rule with that of 

the HDET \W]. HDET has only one scale pp = fJ- and the quark loop contributions 

related to four-quark interaction is of (9(1) so all such loop contributions should be 

#'^There is ambiguity for the overall factor, and thus one may replace ^^2pp/A^" by pp/A^ in 
Eq. 
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summed. However, the same fermionic loop contributions are of 0{{2pp/ K^"^^^) in 
our theory. Though I'ppj is not related to Q directly, it is a small quantity at 
moderate densities which satisfy "pp ^ A^. So Lat is also a good counting parameter 
and one-fermionic loop calculation is enough for moderate density regions, in which 

Q < < Ax (4-5) 

is satisfied. As density increases and pp becomes larger, we should sum the all 
fermion loop with the same v as we should do in HDET. Furthermore, we may have 
to include other degrees of freedom which become light in the high density region. 



4.2 Higher order terms 

One can easily check that all terms listed in the previous section are the leading 
order interactions with z/j = 0. Here we list higher order terms with z/j = 1 for the 
next order calculations. 



Terms with two derivatives 



Possible terms with two derivatives are given by 



(4.6) 



where the covariant derivative acting on is defined in Eq. ()3.1ip . Note that 
the term ^! [7* , 7-'] [Di , Dj\ ^ exists. But this is rewritten as 



7 , T 



(4.7) 



where 



is the field strength defined from the 1-form ay^. We will list this term later. 

*'*Note that the term ^'7*1)^751)°^' is prohibited by the parity invariance. 

#^It should be noted that the free fermion model provides Ci = C2 = C3 = AttF^/h ^ 1. 
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Terms with one meson field and one derivative 

Possible terms including one derivative and one field are given by 



47rF* 



where the derivative acts on only ua field as 



(4.9) 



(4.10) 



Furthermore, as we showed in Eq. ()4.7|) . we have 

hvo 



■ ,T, 

t -. — 



7% 7^ 



47rF* 

TV 

Note that there are identities 



(4.11) 



-Fmi'("V') = i [o^A^l , 01 Au] + , 



(4.12) 
(4.13) 



where A^j^u and V^jy are defined as 



A, 



1 r 

2 
1 

2 



7^ — £ 



1 r 



(4.14) 



with TZ^u and £^,y being the field strengths of the external gauge fields TZ^ and C^. 
So the hvo-teivd in Eq. 1)4.111) can be expressed by the term 



47rF^ 

and the last term in Eq. ()4.16p below. 



(4.15) 
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Terms with two meson fields 

Possible terms with including two of a a field are given by 

V, + ^^aW^ + ^ + , «M ^ • (4.16) 



Six Fermi interaction 

Simplest terms expressing the six Fermi interaction are easily obtained by multi- 
plying each term in Eq. dSH by ^^^) or ^V^#7/.75^) like ^^^)'. 
Other possible terms are expressed as 

le^^o^nV^ |^fr^^75^| |^^'7"75^| l^f'7'375^ 



4vr(F*) 



+ I E ^T'^Tsrvf I I Y: ^7°75r^ ] | E ^7^5^ ] • (4.17) 



Vp I \ Vp I \ Vp 



5 Vector current correlator 

We begin to study our EFT with vector current correlation function. In the present 
case, the vector correlator is simply expressed by the two-point function of the vector 
external field defined from the external fields 7^^ and introduced in Eq. ()3.5|) 
as 

v, = ^{n, + c,) . (5.1) 

Applying our power counting ()4.4p to V^-Vy two-point function denoted by Uy" 
where En = and Ee = 2, we have 

z/ = 2(L-L^) + 5](^rf, + y + e,-2) . (5.2) 
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Figure 1: Diagrams contributing to the Vi_rVi' two-point function: (a) Fermionic 
one- loop contribution with v = Q and Nl — 1; (b) Tree contribution with v = 2 and 
Nl = 0; (c) Pion one-loop contribution with v = 2 and Nl — 0. Here denotes the 
vertex with Vp and • the momentum dependent vertex. 



So the leading loop correction is the fermionic one loop contribution from the 
covariant kinetic term ()3.1U|) . i.e., u = 0,L]\r = 1 (see Fig. Ufa)): *^ 

dH 1 
2(27r)4 I-Vf ■ {I - ViP) - ^e(/o - ViPo)] 

(5.3) 



Vp 



[-Vp ■ {I + r]2P) - ^e(/o + ri2Po)] 

where = (1, 0) and rji and 772 are constants satisfying 771 + ^72 = 1- In this form 

the integration over Zq is finite. Performing the integral over Iq, we obtain 

#6There are no contributions with 1/ = Ln = 0. The contributions with 1^ — 2 and Ljv = are 
expressed by the sum of the tree diagram and the pion one-loop diagram. The tree diagram given 
in Fig. n^b) is constructed from the Lagrangian of the following form: 



where V^''' is defined in Eq. H4.14II . The pion one- loop diagram in Fig.QJc) is constructed from the 
leading order Lagrangian in Eq. (|3.7(l . 
^^In the present integration, we cannot make the shift of the integration momentum in the 

direction of the Fermi velocity. If we were able to make it, we would be able to show that the 
integration in Eq. I|5.3|l vanishes. Then, one might worry that the result of the integral depends on 
the routing of the momentum in each fermion line, which would imply that the final result depend 
on rji and 772. However, as we will show below, the dependence on 771 and 772 will disappear in the 
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Vp 

X / ^ / |j {-v^ ■ {I- mp}) - {-v^ ■ {1 + V2P})] , (5.4) 



Vf ■p + iepo 



where 9{x) represents the step function defined by 



/ X 1 1 for X > , , , 

9{x) = { ' (5.5) 

for a; < . 



In Eq. ()5.4|) we replaced the integration over the spatial momentum as 

d^f r cPl, r dk 



(27r)3 J (27r)2 7 27r 



where 

Vp ■ I 



h = 



Vf 



h 



(5.6) 



l^ = l-vp^, (5.7) 
Vf 



with Vf = \vf\- The integration over leads to 

ns;""'(Po,P) = - E yp^p 1 1^7^ ,, '"'l ■ (5.8) 



We here consider the integration over Z^. Since the momentum l±_ lies on the Fermi 
surface, the combination of the integration over l± and the summation over the Fermi 
velocity vf implies that the contribution from overall Fermi surface is included. 
Thus, it is natural to make the following replacement in the above integral: 



final result, which indicates that the routing of the loop momentum can be done freely at least at 
one-loop level. 
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Then, we have 

nlJ'-lPo.P) = / ^v^v!; ^ '"--l . (5.10) 

It should be noticed that the above contribution is covariant under the spatial 0(3) 
rotation after the angle integration is done. Then, we can decompose this as 

+ prn^(po,p) + prn^(po,p) , (5.11) 

where Pj^ and Pj^ are the longitudinal and transverse polarization tensors defined 
by 

pr = - [g"' - ^) - , 

Pf = 9^{s,,-^^g^, (5.12) 

with p = IpI- By performing the angle integration in Eq. ()5.1()j) . four components of 
Ily'''^'' are expressed as 

KiPo,p) = o, 
Po,p) = — , 

^v{Po,p) = -^^^^Y{Po,p) , 

T-rT/ Pf^F . Pf pI-^fP^v( -\ 



where 



Y{Po,p) 



dilijp vf ■ p 
An Vf ■ p + iepo 

-1 + ^ln^^^±M±!!^ . (5.14) 
2vfP Po - vfP + iepo 
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When the external vector current is electromagnetic, the vector correlator ()5.10|) 
is the photon self-energy and represents the screening effects. Using the formula for 
the static electric potential energy between two static charges qi and q2, 

(27r)'^ + ily[pQ — > 0, p^) 
in Coulomb gauge, the Debye screening mass squared is 

ml = e'^. (5.16) 

which gives the Debye screening mass about 50 MeV at normal nuclear matter den- 
sity. 

The vector correlator Ily'' in Eq. ()5.11|) does not satisfy the current conservation, 
p^Jiy^^" 7^ 0. Like in HDET [131111, this is cured by adding a counter term 

pI [ dn^p ( v'pv^p\ . ^ 



We note that inclusion of the above counter term does not change the Debey screen- 
ing mass in Eq. (j5.16|) . As is well known, the Debey mass comes from the longitudinal 
component Ily, which cannot be expressed by a simple local term in the Lagrangian. 
The above counter term gives a contribution to the spatial component Ily. 



6 Pion decay constants and pion velocity 

In this section we study the correction to the pion decay constants and pion velocity. 
As shown in, e.g., Refs. [1312011211, these physical quantities are extracted from the 
two-point function of the axial vector external field A^, of which Feynman diagram 
is identical to the one in Fig. 1(a) with the external axial vector fields instead of 
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6.1 Fermionic one-loop correction 



Now, let us consider the fermionic one-loop correction to the A^-Ai, two-point func- 
tion. The power counting for the axial vector correlator is the same as the vector 
correlator case (|5.2|) . Then, in the following, we study the fermionic one- loop cor- 
rection to the Af^-Au two-point function. From the interaction term in Eq. ()3.12p . 
the correction with u = and Nl = 1 is expressed as 



5,,n«^'^(po,p) = -trTO]5:2 



K^QUpU a -t- yj _2 



dH 



1 



i{2^Y [-Vf ■ {I - r],p) - te{lo - ViPo)] 
1 



-Vf ■ (/ + ri2P) - ^e(/o + V2P0)] 

(Elr, 



Jab 



pI 



27r'^VF 



i 3 

2 ^t^^i'^F^F 



ao^f'^ + ^A\\gi gj 



VpVp 



+ i^aqk.a\\ {u^gjVp + u'^gjV 



Vf -p 



Vf-p + iepo 



(6.1) 



The contribution to A^-Av two-point function at leading order with v = Nl = 
is expressed as 



(6.2) 

As is done for the vector current correlator in Eq. ()5.11|) . it is convenient to decom- 
pose the sum of the above contributions as 

n7(po,p) = nJ)'^^(po,p) + n«^^(po,p) 

(6.3) 
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From Eqs. ()6.ip and ()6.2|) . two components 11^ and 11^ are obtained as 



(6.4) 



+ I UaoKaW j^-^ALj y{P0,P) 



(6.5) 



where the function Y{pQ,p) is defined in Eq. fl5.14p . 

Using the decomposition of U'^ in Eq. ()6.3p . we can construct the axial vector 
current correlator as 

G7(po,p1 = PrG^APo,P) + PfGlipo,p) , (6.6) 



where 



G'^iPo^p) 
Ga{Po,P) 



A ; 



(6.7) 



We define the on-shell of the pion from the pole of the longitudinal component 
G\ of axial vector current correlator. Since 11^ and 11^ have imaginary parts, 
we determine the pion energy E from the on-shell of the real part by solving the 
following dispersion formula: 







pi Ren^(po, p) - r Ren^(po, P) 



The pion velocity is defined by the relation 

Vt,{p) = E/p 



po=E 



(6. 



(6.9) 



and is then obtained by solving 



Ren^(po = E,p) 
ReU\{po = E,p) 



(6.10) 
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From the expressions of n^(po5P) and Il\{pQ,p) in Eqs. ()6.4p and ()6.5p . Ren^(po 
E,p) = Re[/*]^ and Ren^(po = E,p) = Re[/*/^] are expressed as 

i 

2K 



Re[/:]^ = Ren^(po = E,p) = (F*)^ + ^ (4o4 - ^aoKa\\) X (^) , (6.11) 



Re[/:/^] = Ren^(po = E,p) = F^F^ + 

^2 



Pf 



+ 



2 \ 

All "-AX \ 



2K 



(6.12) 



where K- = \ F^/ F^ is the bare pion velocity and the function X{r) is defined by 



Xir) = 



-1 + —In 
2r 



1 + r 



1 — r 



(6.13) 



Note that we have replaced f,r(p) with \4- in the one- loop contribution in Eqs. (j6.11|) 
and ()6.12|) since the difference is of higher order. 

The damping factor of the pion, 7(p), can be defined from the imaginary part of 
the denominator of the longitudinal axial vector current correlator as 

p 



lip) 



2ReU\{po = E,p) 



Im 



U\{po = E,p)-U:,{po = E,p) 



(6.14) 



By substituting the expressions of Il\{po,p) and Il\{po,p) in Eqs. ()6.4|) and ()6.5p . 
this is rewritten as 



p 2Ren*^(po = ^,P) 27rt;|, 



2 —2 



f^A0K.A\\ 



'^A|| 

4 , 



— K 



A± 



(6.15) 



6.2 Determination of parameters 

We should stress that the parameters of the EFT must be determined from QCD. 
One way to determine them is the Wilsonian matching 1231201; in which the current 
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correlators in the EFT are matched with those obtained in the operator product 
expansion (OPE) around 1 GeV. As is stressed in Refs. j^Hl QB]; when we make 
the matching in the presence of hot and/or dense matter, the parameters of the 
EFT generally depend on the temperature and/or density (the intrinsic temperature 
and/or density dependence). In the present case, however, the cutoff of our EFT is 
around a few hundred MeV at most, so that we cannot make the matching between 
our EFT and QCD directly. So, here we consider the matching between our EFT 
and the HBChPT to determine the parameters included in our EFT. We require 
that the temporal and spatial pion decay constants calculated in our EFT agree 
with those calculated in the HBChPT at a matching point, say pu- 

In the HBChPT at one-loop, the temporal and spatial pion decay constants are 
expressed as jS] 



f{HB)t 

J TT 

f(HB)s 

J TT 



1 + 



2p 



1-^ 

J IT 



C2 + C3 



C2 - C3 



a: 
a: 



(6.16) 
(6.17) 



where f-,, is the pion decay constant at zero density. 

The matching conditions are given by equating Eq. ()6.11|) with Eq. ()6.16p and 
Eq. (I?n^ with Eq. (IHTfj) at p = Pm'- 



P=PM 



(HB)t 



P=Pm 



J TT 



(HB)t f (HB)i 



J TT J TT 



P=PM 
P=PM 



(6.18) 



The above conditions are not enough to determine all the parameters, so we adopt 
the following ansatz to determine the values of kao, i^a\\ and ka± [see footnote |#1| : 

rriN 



-9A 



(6.19) 



where qa is the axial coupling. Here we use the following values of m^r and qa in 
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the matter free vacuum: *^ 



m^ = 939MeV, c/a = 1-267 



(6.20) 



In the approximation ()6.19jl . the temporal and spatial pion decay constants in 
Eqs. (j?rTT|l and (jFTT^ become 



^2 



n2 ^2 



2t;^o P 



(6.21) 
(6.22) 



It should be noted that the damping factor of pion vanishes when we take the ansatz 

(6.23) 



in Eq. (jfTTIH) with /i = Jfp + 



7(p) = . 



Now, by using the matching conditions in Eq. ()6.18p together with the ansatz in 
Eq. ()6.19j) . we determine the values of the parameters and for given matching 
density pu- We take the values for the low-energy constants in the HBChPT as 
C2 = 3.2 ± 0.25 GeV"^ ^ EH and C3 = -4.70 ±1.16 GeV"^ 021 E^. Substituting 
these values into Eqs. ()6.16|) and ()6.17|) . we can reduce them to 



.{HB)t _ . 
f(HB)s _ r 



1 - —(0.26 ±0.04) 
Po 

1 - —(1.23 ±0.07) 
Po 



(6.24) 
(6.25) 



where we take the pion decay constant in the vacuum as = 92.4 MeV and the 

nucleon mass as Eq. ()6.2()|1 . and po denotes the normal nuclear density. When we 

"^^In the high density region, niM should be replaced with to^. Here we consider the density 
region up until around the normal nuclear matter density, so that the vacuum mass will give a 
good approximation. The experimental date shows that there is a 10 % decrease in /^r due to the 
in- medium modification J7|. Thus if decreases by 10 %, our results may be changed by a few 
%. 
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take the matching density as Pm/ = 0.2, 0.3 and 0.4, the values of and are 
obtained by solving the matching conditions Eq. ()6.18|) as follows: 



6.3 Density Dependence of the Pion Decay Constants and 
Pion Velocity 

In this subsection, we study the pion decay constants and pion velocity in dense 
matter. Strictly speaking, the parameters of the EFT generally have the intrinsic 
density dependence. In our construction of the present EFT, we expand the La- 
grangian by and refer to the coefficient of each power of 1/yU as the parameter. 
Then, the density dependences of the parameters are moderate. As a result, the 
density dependences of the physical quantities are dominated by the dense effect 
from the fluctuating nucleon loop in the density region not so much higher than 
the matching point pm ~ (0.2 — 0.4)po- Based on this, we use the values of the 
parameters obtained in the previous subsection to obtain the density dependences 
of the pion decay constants and pion velocity through the dense loop correction 
from the fluctuation mode. We show the density dependences of the temporal and 
spatial pion decay constants in Fig. |21 Until the matching density pm, the behavior 
of physical quantities is described by the HBChPT. The HBChPT breaks down at 
the density where vanishes, around O.Spo- Before the HBChPT breaks down, 
we should switch the theory from the HBChPT to our EFT in which the dominant 
density dependences of the physical quantities come from the nucleon fluctuating 
near the Fermi surface. 

As one can see from the expressions of Eqs. and (j6.17|) . the density de- 

pendence of the pion decay constants in the HBChPT is proportional to p'^p. This 




(88.3 MeV, 69.6 Me V) 
(87.1 MeV, 59.1 MeV) 
(87.3 MeV, 50.6 MeV) 



for pa//po = 0.2, 
for pa//Po = 0.3, 
for pm/Po = 0.4. 



(6.26) 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



0.2 0.4 0.6 0.8 1 1.2 1.4 



(a) (b) (c) 

Figure 2: Density dependences of the pion decay constants with the matching point 
Pm/ Po — 0.2 (a), 0.3 (b) and 0.4 (c). The horizontal axis shows the value of p/poi 
and the vertical axis the value of the pion decay constants normalized by the vacuum 
value. The vertical solid line shows the position of the matching density, Pm/po- Iii 
the low-density region, p < pM, the thick solid lines denote fl/ f-n (upper line) and 
f^l J-n (lower one) obtained from the HBChPT, and the dashed lines denote fH f-,i 
(upper line) and f^/f-w (lower one) obtained from our EFT. In the higher density 
region, p > pM, we used thick solid lines for the predictions obtained from our EFT 
and dashed lines for those from the HBChPT. 

is the reflection that the nucleoli inside Fermi sphere contributes to physical quanti- 
ties. While the pion decay constants in our EFT have the density dependence of p'j^ 
see Eqs. (KT^ and ^TTBf ]. *^ which is a consequence that our EFT includes the 
effects of the nucleon near Fermi surface. Thus the density dependence of physical 
quantities becomes gentle and then the available region extends as shown in Fig. |21 

As well as the HBChPT, our EFT also breaks down at the density where 
becomes zero. Figure |21^c) shows that the EFT is not applicable at p ~ 1.3po- Even 
if is still finite, the EFT is not valid at the density where becomes larger than 

because of the inconsistency with causality. From Fig. |21^a) and (b), we find that 
the EFT breaks down at p ~ 1.3po and p ~ 1.8po, respectively. 

Here in order to check the availability of our EFT, we compare one-loop cor- 
rection with tree contribution. Looking for the density where the ratio of one-loop 

'^^The baryon chemical potential p and function X{vf /V-,^) have the implicit pp dependences. 
However these density dependences are small as compared with p^p. 
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correction to tree contribution becomes 0.5, we find such density as follows ^^^i 
p ^ 1.2po (for Pm/pq = 0.3) and p ^ 0.8po (for Pm/po = 0.4). As we mentioned 
above, Fig. El^a) shows that our EFT breaks down at p ~ 1.3po. Then in the case 
(a) our EFT becomes invalid before perturbative expansion becomes unreliable [The 
ratio at p ~ l.Spo in this case is 0.2 — 0.3, which is still small]. Thus in any case our 
EFT describes physical quantities until around the normal nuclear density po. 

Note that we just calculate u = and Ln = 1. We need to sum all the fermion 
loops with the same u as the density increases. 1.3po indicates how far we can keep 
Ljv as a counting parameter, not how far we can go with counting parameter u. By 
summing fermion loops with 4- Fermi interactions, we can formally extend our theory 
much further. In the high density region, however, the intrinsic density dependence 
of the parameters of the EFT will be important. It is crucial for studying the 
density dependences of the physical quantities that we include the intrinsic density 
dependences of the parameters in some way. Furthermore, we may have to include 
other degrees of freedom which becomes light in the high density region. 

From Eq. ()6.10|) and the above results on /* and we can determine the pion 
velocity Vt^- Figure El shows the density dependence of v.,,. The pion velocity becomes 
smaller with increasing density in low density region. Since the low-energy constant 
C2 in the HBChPT is comparable with tree contribution, the pion velocity largely 
decreases. In the region described by our EFT, the pion velocity scarcely change 
with respect to density. This behavior is quite similar to the one predicted by the 
skyrmion approach in the moderate density region |33]. 
*-'^°The ratio for fl is almost same as that for 
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(a) 



(b) 




Figure 3: Density dependence of the pion velocity with the matching point pul Pq — 
0.2 (a), 0.3 (b) and 0.4 (c). The horizontal axis shows the value of p/po, and the 
vertical axis the value of the pion velocity. The position of the matching density, 
PmIpqi is indicated by the vertical solid lines. For p < pM, the thick solid line 
denotes the prediction of the HBChPT and the dashed line the one of our EFT. For 
p > Pm, on the other hand, we used thick solid line for our EFT and dashed line for 



the HBChPT. 

7 Mass (-like) terms and condensate 
7.1 Two mass terms 

In section 2, we have deleted the original mass term. But the transformation prop- 
erty of \E' given in Eq. ()3.9|) does not prohibit the existence of the term Then, 
one might think that we should include the mass term of the fluctuation field \E' as 



where fi denotes the mass parameter. Since ^"if = \E'^\E', the effect of ft should be 
included into the definition of the chemical potential fi and this term should be 
omitted in the Lagrangian. 

Furthermore, we have another type of mass term which comes from the existence 
of the explicit chiral symmetry breaking due to the current quark mass. In the ordi- 
nary chiral perturbation theory the current quark mass is included by the vacuum 
expectation value of the scalar external source field S. For two flavor case this is 



(7.1) 
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given by 

(5) = A< = f 1 . (7.2) 

In the present analysis we work in the chiral hmit, so that we take 

{S) = 0. (7.3) 

This scalar source field S combined with the pseudoscalar source field V has the 
following transformation property under the chiral symmetry: 

S + zV-^gL{S + zV)g^j, , (7.4) 

where gi and g^ are elements of chiral SU(2)l,r groups. To include S + iV into the 
Lagrangian, we define 

M=iHS + iV) . (7.5) 

This TVf transforms as 

M h{n, g^, gL)-M- h\TX, g^, c/l) • (7.6) 

Then, the term including this A4 field is written as 

Cm = T.Cm^[M+M^]^ , (7.7) 

Vp 

where Cm is a dimensionless real constant. Note that, for constructing the above 
term, we used the parity invariance after summing over the Fermi velocity vp (and 
integrating over x as usual): The transformation properties of and S under parity 
are given by 

"^^^^We define the chemical potential /i at the chiral limit, so that we include the term in Eq. ()7.7|1 
into the Lagrangian. 
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M{xo,x) M\xo,-x) . (7.8) 

Then, 

J d'^x'^{xo,x;vF)M.{xo,x)^{xo,x;vF) 

Vp 

^ y d'^X^{xo^—X]—VF)M^{xQ,—x)'^{xQ,—X]—VF) 

Vp 

= ^J d'^x'^{xo,x;vF)M''{xo,x)'^{xo,x;vF) ■ (7.9) 

Vp 

Since the current quark masses are assigned as 0{Q'^) in the ordinary chiral per- 
turbation theory, it is natural to assign 0{Q^) to the field A4. Then, we count the 
term in Eqs. jlTj) as 0{Q^): 

Cm ~ OiQ') , (7.10) 
or equivalently this term carries z/j = 3 in the power counting given in Eq. ()4.4|) . 



7.2 Quark condensate 

The quark condensate can be read from the one-point function of the external source 
field S. In the present paper we decompose S as 

^ = ^^'+ E ^^T,^ (7.11) 
iv f 1 

/ a=l 

where Ta is the generator of SU(A''/) normalized as tr[TaTfe] = ^6ab- Then, the 
one-point function of provides the quark condensate per one fiavor. From the 
interaction in Eq. ()7.7|) . one-loop correction to one-point function of is expressed 
as 
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A 



dH 



(7.12) 



This integral is divergent, so that we here define the regularization to perform the 
integral We first replace the momentum integral as in Eq. ()5.6|) . Then, we introduce 
the cutoff for Z|| and Iq as — A|| < /|| < A|| and — Aq < Iq < Aq. By using this 
regularization method, Eq. ()7.12j) is written as 



Vf 



Performing the Iq integral, we obtain 



Vp 



X 



-e{-VFi\Oe{A^-\vFi\\\) + 



+ ^ {in (Ao + {vfI\^{1 - te)) - In (Aq - {vfI{){1 - ^e))} 



(7.14) 



By taking the ^ oc limit, this expression is reduced to the following finite form: 



_ ^ /■ (fL r dill ^, . 
^^-^^-5/(2^/2^^(^1 



- l/ii 



-vfI\\) + - 



(7.15) 



Now, performing the /ii integral, we obtain 



A = . 



(7.16) 



This implies that the number density of the fiuctuation mode at one-loop level is 
zero in chiral limit. It is a kind of a symptom that we choose a proper vacuum in 
matter. 
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8 Discussions 



As the density of the nuclear medium increases, the Fermi momentum pp emerges 
as an additional scale of the system, which deserves to be treated differently from 
the small scale Q of HBChPT, for a more transparent and better chiral expansion 
scheme. In this work, we have considered the cases where is much larger than Q 
but still much smaller than the chiral symmetry breaking scale, Q <^ pp <^ A^- 

We studied the vector and the axial vector correlators in the nuclear medium, 
and we could calculate the Debyc screening mass, pion velocity and the modification 
of the pion decay constant. The density-dependence of those quantities is found to 
be highly different from that of the HBChPT, even at the moderate density region. 
Instead of having the leading linear density dependence as was the case in HBChPT, 
we have observed rather mild dependence. This may be understood by a following 
simple consideration: At very low density, the response of the external probe would 
be proportional to the number of particles in the system, or to the volume of the 
Fermi sea. When density increases, however, excitations are limited around the 
Fermi surface due to the Pauli blocking, and the response would now be proportional 
to the area of the Fermi surface. 

The present study enables us to describe the nuclear matter up to the density 
comparable to po, the normal nuclear matter density. 

To extend the applicable region of the theory, we need to do RPA type of cal- 
culations to sum all the contributions with respect to pp/A-^, for a given order of 
Q/A^. This step will result in a very rapid convergence for the processes where 
Q/A^ -C 1, regardless of the density. The actual calculation can in fact be done 
rather easily. But this step makes our results to be dependent on the coefficients 
of the four-fermion interactions. And thus we need a more complete study includ- 
ing the renormalization of those coefficients to have any meaningful predictions. In 
addition to this, as the density increases, the intrinsic density dependence of our 
theory can play a crucial role, and thus should be properly treated. The notion of 
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the decimation [33], and the renormahzation-group running and additional matching 
conditions might be necessary to address this issue. Studies on these subjects are 
left for the future work. 

In the high density region in which the RPA is needed, the perturbative structure 
of the chiral effective theory will be most likely modified with the necessity of choice 
of the relevant degrees of freedom that we play with. Indeed, we may need to 
introduce the light vector meson degrees of freedom as claimed to do for the matching 
with the QCD in the phase of the Vector Manifestation ^H]- Even with the chiral 
partner of the pion is incorporated, the perturbative scheme should stay unchanged 
as shown in the counting rule. 

However, when the chiral phase transition is approached, we have to reconsider 
the change of the representative scales of the medium from what is considered in this 
work. Namely, the 47r/^ is reduced eventually to be zero, and thus it cannot serve 
as the "heavy scale" of the system. In this case, the relevant expansion parameter 
would be the ratio of the typical momentum transfer compared to the chemical 
potential, Q/yU, which results in a theory with the same scale dependence as the 
HDET [16^. The transmutation of the scales could be realized by considering the 
second decimation as is suggested by Brown et al [HB] . 
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